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X-ray flares reveal mass and angular momentum of the Galactic 

Center black hole 
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Abstract. We have analysed the light curve of the two brightest X-ray flares from the Galactic Center black 
hole, one flare observed by XMM-Newton on October 3, 2002 (Porquet et al. 120031 . and the other flare observed 
by Chandra on October 26, 2000 (Baganoff et al. 1271011 . The power density spectra show five distinct peaks at 
periods of ~ 100 s, 219 s, 700 s, 1150 s, and 2250 s common to both observations within their estimated measurement 
uncertainties. The power density spectrum of the recently reported infrared flare of June 16, 2003 (Genzel et al. 
2003) shows distinct peaks at two, if not three, periods (including the 1008±120s infrared period), which are 
consistent with the X-ray periods. The remaining two periods could not be covered by the infrared measurements. 
Each period can be identified with one of the characteristic gravitational cyclic modes associated with accretion 
disks, i.e. either Lense-Thirring precession, Kepler orbital motion and the vertical and radial epicyclic oscillation 
modes, in such a way that a consistent value for the black hole mass of M BH = 2.72±° i 2 x 1O 6 M and the 
angular momentum a = 0.99391q'qo74 ^ s obtained. The available data on Mbh derived from studies of the orbital 
motion of the S2 (SO-2) star (Schodel et al. l2T)CJ2l Ghez et al. 2003) agree with our findings. Finally we discuss 
some implications of the fairly high value for the angular momentum. 

Key words. Galaxy: center - X-rays: individuals: Sgr A* - X-rays: general - Radiation mechanisms: general 



1. Introduction 

The region towards the Galactic Center (GC) has been 
recently resolved in X-rays by Chandra measurements 
(Muno et al. I2U031 Baganoff et al. I2UU5|) . A source was 
identified as X-ray counterpart of the GC black hole 
SgrA*. The quiescent X-ray luminosity turns out to be 
extraordinary low (2.2x 10 33 ergs -1 , Baganoff et al. 2003) 
for a mass accreting black hole, given its mass which has 
been measured through the motion of the SO-2 star or- 
biting the black hole. Mass values range from 2.2xlO 6 M 
(lower limit of Schodel et al. l2UU2)l . 3.6 x 10 6 M Q (Genzel 
et al. H003I> to 4.1±0.6 xlO 6 M for a distance of 8 kpc 
(Ghez et al. [2003 ) . Various models to explain the low lu- 
minosity have been proposed and discussed quite contro- 
versially including low accretion rates and radiatively in- 
efficient accretion flows, for example. For a brief overview 
see for instance Melia & Falcke {20010. Baganoff et al. 
PI) . Porquet et al. ({2731)5)) and Yuan et al. (|27JTJ5j) . The 
quiescent state is interrupted by occasional, sometimes 
very bright, flares. The first bright flare with a peak lu- 
minosity of lx 10 35 ergs _1 (Baganoff et al. 2001) was ob- 
served with the Chandra ACIS-I instrument on October 
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26, 2000. An even brighter X-ray flare with a peak lu- 
minosity of 3.6xl0 35 ergs -1 occurred on October 3, 2002 
and was caught with the XMM-Newton EPIC instruments 
(Porquet et al. 2003 ) . ft is generally accepted that the ra- 
diation is from or at least channeled through the accretion 
disk orbiting the central black hole. Some models for the 
creation of flares have been brought forward, e.g., Markoff 
et al. I|2001[) suggest extra electron heating near the black 
hole, Liu & Melia ( 2002) propose a sudden enhancement 
of accretion and Nayakshin & Sunyacv ( 2003 ) prefer stars 
passing the disk. We think that the high luminosity in- 
volved in the flares suggests that somehow the inner re- 
gions of the accretion disk play a major role, basically 
because of the increased gravitational potential. Sporadic 
events like flares are then likely to take us closer to the 
black hole and enable us to discover more easily periodic 
or quasi-periodic changes in the light curve of a flare, like 
the Kepler frequency of the innermost marginally stable 
orbit, and others, by which mass and angular momentum 
of the central black hole could be determined. That this 
idea might have some bearing has been recently demon- 
strated by the discovery of a 16.8±2 min period in two 
infrared flares (Genzel et al. 120031 . Genzel et al. tenta- 
tively assign the period to the length of the marginally 
stable orbit and, adopting a mass of 3.6 x 10 6 M Q , they 
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calculate an angular momentum of a = 0.5 for the black 
hole. 

In section 2 we show the light curves of the two bright- 
est X-flares, which we have analysed including but sepa- 
rated from the time sections before and after the flare 
proper, as well as two very long (11 hours and 46 hours) X- 
ray light curves representing the so-called quiescent state. 
Like the two flares one event was observed by XMM- 
Newton and the other event was covered by Chandra. 
These two observations were primarily analysed for cross- 
checking with the flare characteristics. Section 2 also con- 
tains the power density spectra (PDS) of the relevant light 
curves. Five groups of characteristic frequencies with at 
least two detections per group have been found. In sec- 
tion 3 we describe the relation of these periods with the 
cyclic gravitational modes thought to be possibly created 
in accretion disks, i.e. the Lense-Thirring precession fre- 
quency, the Kepler frequency and the vertical and radial 
epicyclic frequencies. In section 3 we also describe how 
we derive the values for Mbh and the angular momentum 
a (the Kerr parameter) from the frequencies measured. 
We include in our analysis the power density spectra of 
the two infrared flares published by Genzel et al. |2003 ) . 
Finally we discuss a few implications which result from 
the fairly high value of a we have obtained. 



2. Observations and analysis 

2.1. Light curves 

Figure^ shows the light curve of the October 3, 2002 flare 
observed by XMM-Newton and published by Porquet et 
al. (2003) with a slightly different time binning. Figure^) 
shows the light curve of an ~11 hours quiescent period ob- 
tained by XMM-Newton on February 26, 2002. Figure 
displays the light curve of the very first flare reported from 
SgrA* dated October 26, 2000, observed with Chandra 
and published by Baganoff et al. (201) 111. We choose a 
slightly different energy band and a different binning. 
SgrA* was observed for 46.5 hours on May 5/6, 2002 with- 
out a bright flare occurring, and the Chandra light curve 
is shown in Figure^. Because of the numerous peaks the 
source is apparently not in a true quiescent state but ex- 
hibits quite a number of smaller flux increases or small 
flares. The essential observational details are summarized 
in Tabled The XMM-Newton data have been proprietary 
PI data, which by now are publicly available. The Chandra 
data have been extracted from the public Chandra archive 
(level 2 processed event list, provided by the Chandra X- 
ray Center). We restrict our analyses to just these four 
data sets because they cover the two brightest flares and 
the other two observations have the longest exposure of 
the so-called quiescent state. For the timing analysis we 
added PN, MOS1 and MOS2 data in one set using counts, 
which were for XMM-Newton extracted from a circle of 
10" radius with energies between 2.6 - 10 keV. For the 
Chandra data the extraction radius was 1.5" around the 



Table 1. Summary of observations analysed (c.f. 
Figs.UEJ). Instr. I is EPIC PN+MOS1+MOS2 on XMM- 
Newton, and instr. II is ACIS-I on Chandra. Date is in 
day/month/year. 



Instr. 


observ. 


total 


flare 


preced. 


sampl. 






time 


time 


time 


time 




date 


(min) 


(min) 


(min) 


00 


I 


03/10/02 


254.1 


56.33 


197.8 


2.6 


I 


26/02/02 


667.2 






2.6 


II 


26/10/00 


589.3 


115.3 


372.4 


3.241 


II 


25/05/02 


2804.2 






3.141 



archived source position, and the energy band is 2.0 - 8 
keV. 

Because of the relatively broad point spread function of 
XMM-Newton the extraction radius of 10" contains emis- 
sion from more objects than just SgrA*. The Chandra 
image reveals another four or five additional sources as 
well as diffuse emission. Baganoff et al. (2003; estimate a 
net count rate of 5.74±0.40 x 10~ 3 counts -1 for SgrA* 
in its quiescent state. The power law spectrum has a pho- 
ton index of T q = 2.7^J'g and an absorbing column density 
Nn, q =9.8^3qX 10 22 cm -2 . Using the best fit values we es- 
timate a SgrA* quiescent state count rate of 0.029 ± 0.003 
counts" 1 for MOS 1+MOS2+PN, so that the conversion 
factor is 5 ± 0.5 between the two instrument combinations. 
The mean EPIC count rate of Sgr A* in its quiescent state 
amounts to 0.16 counts" 1 (c.f. Fig.^J)), which means that 
the background sources are likely to contribute on average 
0.013 count s" 1 or 82% of the total during the Sgr A* qui- 
escent state, provided that neither the background sources 
nor SgrA* have changed significantly compared to the 
October 26, 2000 observation. When we subtract the count 
rate of the background from the flare observation (fig.[IJi) 
the peak count rate during the time preceding the actual 
flare has a level of 0.125±0.03 count s , which means that 
the precursor section is ^4 times brighter than SgrA* in 
its quiescent state. It is interesting that flares observed in 
the infrared by Genzel et al. (|2003|) show about the same 
ratio of the flare peak intensity and the quiescent level. 
This suggests that the X-ray precursor and the infrared 
flares originated in the same process and that, given the 
huge X-ray flares, much brighter infrared flares still await 
their detections. 

Conversion of instrument count rates to flux density 
requires knowledge of the source spectrum and column 
density of the photoelectrical absorption. This is available 
for the quiescent state measured by Chandra (Baganoff 
et al. E133, the October 26, 2000 flare measured by 
Chandra (Baganoff et al. USD and the October 3, 2002 
flare measured by XMM-Newton (Porquet et al. 12003(1 . 
Since both the spectrum and the absorption column den- 
sity differ among the three observations we provide con- 
version factors for the EPIC count rate and the ACIS-I 
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Fig. 1. EPIC light curves (MOS 1+MOS 2+PN) of the 
XMM-Newton observation of October 3, 2002 (upper 
panel, Fig. la) and the February 26, 2002 observation 
(lower three panels, Fig. lb). Error bars indicate la un- 
certainties. The horizontal line in the lower three panels 
corresponds to the mean count rate level. Arrows mark 
peaks associated with a 2178 s periodic signal. 

count rate for each of the three states. These are: 9.5 ^Jy 
(EPIC counts/s)" 1 , 47.6 ^Jy (ACIS-I counts/s)- 1 (quies- 
cent); 0.5±J:| fiJy (EPIC counts/s)" 1 , 2.8±|;J n.Jy (ACIS-I 
counts/s)" 1 (Chandra flare); 12.3 /iJy (EPIC counts/s) -1 , 
71.8 ^Jy (ACIS-I counts/s)" 1 (XMM-Newton flare). The 
relatively large uncertainty of the conversion factor for the 
Chandra flare are due to the uncertain spectrum and col- 
umn density. The energy flux densities are given for E = 
1 keV. 

2.2. Power spectra 

From the light curves shown in Figuresnand[2 power den- 
sity spectra (PDS) have been created by a Fourier anal- 
ysis of six data sets. Two sets cover the time section of 
just the flare, one each for XMM-Newton (Fig. [3]}) and 
Chandra (Fig. 0Jl) ; another two sets cover the time sec- 
tions before and after the flare (Figs. |3Jd and @J>). For 
these PDS the data of the flare proper have been removed 
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Fig. 2. ACIS-I light curves of the Chandra observation of 
October 26, 2000 (upper panel, Fig. 2a) and the May 25, 
2002 observation (lower three panels, Fig. 2b). Error bars 
indicate la uncertainties. 

and replaced by data with a mean flux identical to that 
of the time preceding the flare assuming a Poissonian sta- 
tistical distribution. The final two sets correspond to the 
two observations of the quiescent level (Figs.|3]and[f)J). We 
define the power spectral density psd ra at a frequency f n 
and wave number n as psd„ = (a^+b^)/(2x At 2 ) x T; a„ 
and b„ are the Fourier coefficients, At is the sampling or 
binning time and T is the total observing time (Table ^) . 
The Fourier coefficients are in units of ACIS-I counts -1 
for the Chandra observations and EPIC counts -1 for the 
XMM-Newton observations. Because of the different effi- 
ciencies the psd n 's are expected to generally differ by a 
factor of ^25±5 for the same luminosity level given the 
energy spectrum described in section 2.1. 

2.3. Frequencies 

In general the six PDS look very much alike. There is a 
high frequency component which is more or less rapidly 
growing with increasing frequency, which reflects the 
noise introduced by the low counting statistics. We note 
that the smallest possible binning has been chosen so 
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Fig. 3. Power density spectra of the October 3, 2002 flare Fig. 4. Power density spectrum of the October 26, 2000 
(upper panel, Fig. 3a) and its precursor section (lower flare (upper panel, Fig. 4a) and its precursor section (lower 



panel, Fig. 3b); XMM-Newton measurement. 



panel, Fig. 4b); Chandra measurement. 



that we deal with just no or one count per bin. At low 
frequencies the PDS are dominated by a few high peaks, 
which summarize the mean shape of the light curves, 
averaged over periods which are just a few times shorter 
than the whole track. Clearly separated from these two 
regimes is a low noise mid-frequency section between 
~0.7mHz and ~7mHz for the flare observations (Fig. [3ji 
and^Ji), and 0.7 mHz and 2.5 mHz for the flare precursor 
sections (Figs. 0h> and EJ)). In this frequency band we 
see two distinct PDS peaks at 1.426 mHz and 4.562 mHz 
(Fig.|2^; labels 3 and 1, respectively; XMM-Newton flare) 
and similarly at 1.445 mHz and at 3.902 mHz (Fig. 0^; 
labels 3 and 1, respectively; Chandra flare), standing 
well above the noise level in the corresponding frequency 
bands. In the Chandra flare we see a third peak at 
10.41 mHz (label 0), which, however, is located at the 
beginning of the climbing high frequency noise section 
(Fig. Hi). The XMM-Newton flare PDS shows a similar 
close-by peak at 9.12 mHz (label 0), although at a fairly 
low power level. Like in the Chandra case this frequency 



is at the high frequency edge of the low-noise band. The 
existence of this pair of frequencies might be doubtful 
but they should be noted basically because they seem to 
appear in two different observations. These are the only 
pairs in the flare section. Absent in the XMM-Newton 
PDS but fairly prominent in the Chandra PDS (Fig. 
is a peak at ~2mHz (label 2). 

The PDS of the data preceding the flare also show dis- 
tinct peaks in the low noise regime at 0.853 mHz (Fig. 0t>, 
label 4; XMM-Newton) and 0.895 mHz (Fig. Hj, label 4; 
Chandra); a further peak appears in the XMM-Newton 
PDS at 0.459 mHz (label 5). The power at this frequency 
is so high that the equivalent period of 2178 s can be iden- 
tified by peaks in the light curve (c.f. arrows in Fig.^,). A 
PDS peak associated with or close to 0.459 mHz is clearly 
absent in the Chandra flare precursor observation but 
there is a peak at 0.434 mHz in the Chandra flare (label 5). 
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Fig. 5. Power density spectrum of the February 26, 2002 
quiescent state; XMM-Newton measurement. 
0.18 
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Fig. 6. Power density spectrum of the May 25, 2002 qui- 
escent state; Chandra measurement. 

The PDS is a function of discrete frequencies given 
by f„ = y with 1< n < int( 2 T At ) and At the binning 
size. Therefore, there is in principle a systematic rela- 
tive frequency uncertainty possible of (4^) n = ±(-^j). 
In Table [3 we summarize the results including frequency 
/, wave number n, period P, power spectral density psd 
and a label ID, which has the same value for frequencies 
close to each other. The same ID's are shown as labels 
in the PDS graphs. We find in the XMM-Newton and 
the Chandra observations five groups of periods, each of 
which is a pair with one member from XMM-Newton and 
one from Chandra. The periods of the members of each 
pair are almost identical, i.e. 110/96s (label 0), 219/256s 
(1), 701/692s (3), 1173/1117s (4), 2178/2307s (5). They 



Table 2. Compilation of outstanding frequencies. ^ The 
psd is given in units of (counts -1 ) 2 Hz" 1 . For XMM 
EPIC count rates and for Chandra ACIS-I, count rates 
are quoted, na: not applicable, i.e., not covered because 
the flare didn't last long enough. The lower section of the 
table contains the frequencies of high psd observed in the 
two infrared flares by Genzel et al. ( 200SJ) . The data have 
been read off their figure 0;. 



State/ 












Instrum. 


f 


n 


P 


psd 


ID 




(mHz) 




(8) 


i) 








flare 








XMM 


9.123 


32 


110 


6.3 





Chandra 


10.41 


72 


96 


0.78 





XMM 


4.562 


16 


219 


15.3 


1 


Chandra 


3.902 


27 


256 


0.78 


1 


XMM 








- 


- 


Chandra 


2.023 


14 


494 


0.46 


2 


XMM 


1.426 


5 


701 


16.4 


3 


Chandra 


1.445 


10 


692 


0.49 


3 


XMM 


na 


na 


na 


na 


na 


Chandra 


0.434 


3 


2307 


1.6 


5 




precursor 






XMM 


0.853 


13 


1173 


1.3 


4 


Chandra 


0.895 


20 


1117 


0.11 


4 


XMM 


0.459 


7 


2178 


1.9 


5 


Chandra 
















IR flare 








IR/16 


4.76 


24 


214 




1 


IR/15 


2.0 


6 


498 




2 


IR/16 


1.39 


7 


733 




3 


IR/15 


1.0 


5 


996 




4 


IR/16 


0.99 


5 


1026 




4 



appear to be not exactly identical but conceding the max- 
imal possible uncertainty they are consistent with each 
other. This coincidence strongly supports their existence 
and suggests that each pair represents the same process. 

Genzel et al. I|2003|l have published the discovery 
of a 16.8±2min period in the two infrared flares ob- 
served on June 15 and June 16, 2003. A look at the 
two published PDS shows that there are more peaks, 
which we have read off from their figure |2t and added 
to Table H (IR/15, IR/16). Except the peak at 321s 
and wave number n = 16, appearing in the June 16 
flare, each frequency has a close-by counterpart in the 
XMM-Newton and/or Chandra PDS adding another 
finding to groups 1, 2, 3 and 4. Of course, the psd at 
214 s in IR/16 is fairly low but it shows up as a local 
maximum. There is no peak at the frequencies of group 
5 in the infrared observations, which however have 
been too short to search for it. Obviously, the infrared 
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frequencies are fully consistent with the X-ray frequencies. 

We note that at none of the frequencies selected from 
the flares or precursors a prominent psd is evident in the 
quiescent level observations. The PDS of the quiescent 
level observations are shown in Fig.[S](XMM-iVewtan) and 
Fig. (Chandra). The XMM-Newton PDS shows a set of 
six prominent PDS peaks between 0.65 mHz and ~2mHz 
centered on — lmHz. This frequency band embraces the 
0.853/0.895mHz peaks of both the XMM-Newton and 
Chandra flare precursors. The psd's at the six frequen- 
cies vary between 1.05 and 1.42 (EPIC count s^ 1 ) 2 Hz -1 
which is close to what is observed as power for the 1173s 
period for the flare precursor (Table |2J . A psd peak at 
the period of 2178 s is not present in the quiescent obser- 
vations. The PDS of the Chandra quiescent level obser- 
vations doesn't show prominent peaks around one mHz 
(Fig.©. But extrapolating the XMM-Newton psd of -1.2 
(EPIC count s -1 ) 2 Hz -1 the equivalent Chandra psd is 
expected to be close to the noise level. 

The lowest frequencies accessible for this investiga- 
tion are covered by the Chandra May 25, 2002 obser- 
vation, with no bright flare though (Figure . Starting 
at —0.1 mHz the psd drops rapidly with frequency reach- 
ing the noise level minimum at —0.4 mHz. The region be- 
tween 0.18 mHz and 0.35 mHz seems to show excess power 
peaking at - 0.24 mHz or 4100 s. In this band the XMM- 
Newton PDS is very similar to the Chandra PDS with 
fairly high peaks at 0.225 mHz and 0.1 mHz, but details 
like in the Chandra observation are not resolved because 
of the shorter exposure (Figure ©. This PSD structure 
is probably not related to the occurrence rate of flares, 
which has been estimated to 1.2±0.6 per day by Baganoff 
(2003), which is a factor of about 17 lower in frequency. 
Excluding the 0.18 mHz - 0.35 mHz band the PDS, with 
20 single psd values remaining, can be fitted by a power 
law with an index of roughly —2.3. We stress that the 
statistics are poor and we don't claim the existence of a 
QPO but the excess between 0.18 mHz and 0.35 mHz looks 
like a QPO structure. QPO's (quasi-periodic oscillations), 
which show up as a band of increased psd above an under- 
lying, much broader PDS, which has the shape of a power 
law, have been observed in a number of galactic binaries, 
containing a compact object like a neutron star or black 
hole. The QPO's are considered to be created in associ- 
ated accretion disks (e.g., van der Klis. I200"n| . As Nowak 
& Lehr l|1998[) and others have pointed out, QPOs aris- 
ing from supermassive black holes in AGN are expected 
to show periods between three hours and years, if a 0.1 s 
period is typical for a galactic black hole and if periods 
scale with mass. Such long periods are very difficult to 
measure for current missions. But for the Galactic Center 
black hole with its fairly low mass QPO periods might 
be shorter. In fact, if we scale the 300 Hz, which is the 
highest QPO frequency observed in the 7 M Q black hole 
microquasar GRO J1655-40 (Remillard et al. 119991 Orosz 
& Bailyn 199*3) we expect a QPO around 0.6 mHz for the 



Galactic Center black hole, which is well within the long 
duration Chandra and XMM-Newton measurements. 

2.4. Statistical significance 

To quantify the significance of a measured psd value we 
have approximated the rather erratic flare light curves by 
some smooth functions which reproduce very well the gen- 
eral shape but do not introduce high frequency variations 
in the regime which contains the frequencies to be inves- 
tigated. We are particularly interested in the frequencies 
of group and group 1. The lower frequencies may be 
too much biased by red noise for which we don't have a 
model. The model light curve of the XMM-Newton flare 
consists of three exponentials, one for the rise, one for the 
top and the third for the decay. The model light curve for 
the Chandra flare is set up by a four piece polygon. Each 
model light curve is normalized to the total number of 
measured counts. The model light curves are binned with 
the same sampling time as the measured data, and the 
model counts per bin are subject to Poissonian statistics. 
For each flare Z = 3x 10 5 light curves have been produced, 
Fourier analysed and the PDS calculated. 

The significance of a power peak is defined via a confi- 
dence level S, which gives the probability that the detected 
power psd nj dotcct is not produced by the noise process, 
i.e. 1 - S = NtrialsXW(psd n > psd n ,dotoct) (van der Klis, 
1989). W is the ratio of the number of powers exceeding 
psdnjdotcct and Z. N tr i a is is the number of frequencies for 
which the PDS has been measured, which is 674 and 1066 
for the XMM-Newton and Chandra flare, respectively. We 
stress that we have not made use of a fixed form for W like 
a x 2 or Gaussian distribution, but we built W by Monte 
Carlo simulations (c.f. van der Klis, 1989). The simula- 
tions show that there are some significant peaks, i.e. for 
group 0: S(96 s) = 0.989 whereas the psd at 110 s of the 
XMM-Newton flare is consistent with a noise created peak; 
for group 1: S(219 s) = 0.993 and S(256 s) = 0.947, which 
means that signals with the latter periods quite likely ex- 
ist. There is just one more pair of close-by frequencies each 
of which shows a significant psd, which are at a period of 
11.31 s with S(11.31 s) = 0.634 in the XMM-Newton flare 
and of 10.68 s with S(10.68 s) = 0.968 in the Chandra 
flare. This very short period should be kept in mind for 
future studies. 

We also tried quite a few different model light curves 
but they do not change significantly the values for S given 
before. 

3. Discussion 

Five characteristic frequency groups have been discovered 
in the X-ray flares with their existence being supported by 
the fact that they have been found in more than one ob- 
servation. We try to identify each of these frequencies with 
one of the oscillation modes which could occur in accretion 
disks surrounding black holes. As Nowak & Lehr ( 1998) for 
instance point out there are four cyclic modes which could 
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give rise to periodic or quasi-periodic changes in accretion 
disks. These modes are based on gravitation but it is an 
open question whether they would create changes of the 
light output at the very same frequencies. However, Nowak 
& Lehr 1|1998[) point out that the mode frequencies pre- 
dominantly depend upon fundamental gravitational fre- 
quencies and are not strongly affected by hydrodynamic 
processes for thin disks. 

3.1. Characteristic accretion disk frequencies 

There are four cyclic gravitational modes associated with 
black hole accretion disks (Nowak & Lehr 1998, equations 
^ to |21 and Merloni et al. 119991 equation 0J, which are the 
Kepler frequency (Ok), the disk perturbation frequencies 
in vertical and radial direction called vertical (fly) and ra- 
dial (fin) epicyclic frequency and the Lense-Thirring pre- 
cession frequency (Olt)- Each frequency depends on the 
central mass M, the angular momentum a and the radial 
distance r from the center. EquationsnjtolHshow the rela- 
tions, for which the standard notation of c=G=l is used. 
Physical length scales are in units of GM/c 2 and angular 
frequencies O are in units of c 3 /GM. r = 1 is defined as 
the gravitational radius r g . 

K = (r 3 / 2 + a)- 1 (1) 

o 2 , = nl (i - 4^ + 2£) ( 2 ) 
o| = o K (l - s + - *£) (3) 

L T = fi K - fi V - (4) 

3.2. Arrangement of frequencies and modes 

Equations H] to 0] show that Ok decreases monotonically 
with increasing r and, on the contrary, for a fixed, both 
fly and R initially rise with r, go through a maximum 
and then decrease with increasing r. fly reaches its max- 
imum at a radial distance denoted by r = rymax and flu 
takes its maximum value at a distance denoted by r — 
r Rmax- For any r and a fixed Ok > fly > R . This means 
that the shortest period is to be assigned to Ok at the 
smallest r possible which is the innermost marginally sta- 
ble orbit at r = r ms defined by R = 0. For each of the 
other frequencies there is a choice between Qy(r = r ms ), 
O v (r = rymax), fly(r = r Rmax ), R (r = r Vmn ) and R (r 
= rRmax). R (r = r ms ) is excluded because R (r = r ms ) 
= per definition. Whether Olt could be generated at 
T > ims is hard to tell, but Olt(?* = r ms ) could approach 
R (V = r ms ) or could even be lower than the frequency 
of any other mode, depending on a. We have chosen three 
radial fix points, i.e. r = r ms , r = r Vroax and r = r Rmax . 
In principle there is no preference for any r except for r 
= r ms but we note that fly and R do not vary a great 
deal with r over a reasonable range of r centered on ry m ax 
or r Rmax , respectively. In the following we use periods, 
i.e. Plt, Pkj Pv and P R instead of angular frequencies 



O. Each observed period is now assigned to one mode at 
either r ms , rymax or r Rmax each resulting in a relation be- 
tween Mbh and a. If these assignments cover all observed 
periods within their measurement uncertainties, a com- 
mon value of Mbh and a should emerge. Group contains 
the shortest period of 110s; if assigned to Px(r = r rms ) 
an upper limit for Mbh < 1.78xl0 6 M Q exists for a — ► 1, 
which is in conflict with the lowest value of 2.2xl0 6 M Q 
published by Schodel et al. 1|2002)) . Therefore this ascrip- 
tion is excluded. The next longer period of 219 s implies an 
upper limit of Mbh < 3.55xl0 6 M Q if assigned to Pk(t = 
r rms ). This is within the range of values for Mbh reported 
in the literature and suggests the assignment to Keplerian 
motion along the marginally stable orbit. 

The only solution of equations to 0] with consistent 
values for Mbh and a is obtained for the following assig- 
nation of period and mode: 
group 1, 219s = P K (r-r ms ) [K]; 
group 3, 692 s = P v (r=r ms ) [V]; 
group 4, 1117s = P R (r=r Rmax ) [R]. 

With this set equation 01 predicts Plt ~ 320 s. Such a pe- 
riod shows up with increased power in the June 16 infrared 
flare (321s, Genzel et al. l2"UU51 c.f. section 2.3). The four 
relations for Mbh versus (1 — a) following from equations^ 
to 01 are drawn in figured which shows an almost perfect, 
single intersection of the four relations. The wide error 
polygon of figure {7\ indicates the maximal range for Mbh 
and (1 — a) without the 321s period and the inner poly- 
gon covers the constraints set by the 321s period. This 
solution does not suggest a specific radius like r ms , ry m ax 
or r Rmax for the remaining period of ~500 s, but it is con- 
sistent with a Kepler period expected close to r=r Rmax . 
The periods of group 5 around 2200-2300 s, cannot arise 
from a region r <r Rmax . However, we point out that this 
period is very close to just twice the period of the radial 
epicyclic mode and might indicate a dominance of the am- 
plitude of every second cycle. We note that the group 
period of ^100 s is close to the first harmonic (m = 2) ex- 
pected for the 'diskoseismology' C-mode (Nowak & Lehr 

MB- 

The best-fit solution to the intersections is Mbh = 
2.72±°;^ x 1O 6 M and a = 0.9939^;^. The errors 
quoted are not la errors but reflect the maximal possible 
uncertainty given by the maximal uncertainties of the pe- 
riod determination. The radii involved are r ms = 1.371 and 
r Rm ax =2.507. The gravitational radius (r=l) amounts to 
4.1xlO n cm. 

Independent measurements of the orbital motion of the 
star SO-2 (S2) around the GC black hole have resulted in 
M BH = 3.7±1.5 x 10 6 M Q (Schodel et al. HUU^i and M BH 
= 4.07±0.62 x 10 6 M Q (R /8kpc) 3 (Ghez et al. 1211)5)) with 
Ro the distance to the black hole. Our result is consistent 
with the Schodel et al. measurement at their 1 a level 
and it also agrees with Ghez et al. 's measurement at their 
2a level and Ro = 8 kpc. Earlier Genzel et al. {1996) 
determined Mbh to 2.5 - 3.0 x 10 6 M Q from measurements 
of dark matter concentration in the central parsec of the 
Milky Way which appears to match our result best. 
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Fig. 7. Relation of Mbh versus (1 — a). Each curve repre- 
sents one period assigned to a specific gravitational mode, 
coded by a label described in the text. The error region 
is enclosed by the dashed polygon. The dotted polygon is 
narrower because of the additional constraints imposed by 
the 32fs period (LT). 

3.3. The high value of a 

The high value of a ~ 1 means that the emission from 
the inner parts of the accretion disk is quite close to the 
black hole. The event horizon is located at rjj = 1.112. 
With r ms = 1.371 the distance to the marginally stable 
orbit is just 3.6 light seconds. Light from matter cross- 
ing the marginally stable orbit is therefore very likely to 
disappear in a fairly short time. Is the PDS peak at ~11 
s the characteristic time scale at which the light is fad- 
ing away? The light curve of the XMM-Newton flare (fig- 
ure^,) shows a deep, very short (< 100 s) cut close to 
maximum light, which might be such an event. A similar 
event is also seen in the Chandra flare (figure ■ The 
arrows plotted in figure point to the maxima of the 
2178 s period, which seem somehow to synchronize with 
the deep cut in the light curve. According to our arrange- 
ment scheme the 2178s period is characteristic for the 
radial epicyclic mode and it looks as if matter is pushed 
across the marginally stable orbit. 

The high value of a « 1 and the proximity of the in- 
ner region to the event horizon also implies that fairly 
low (<1) general relativistic boost factors g have a ma- 
jor impact on the emission because of increasing beaming 
(e.g. Miiller & Camcnzind 2003J. The radiation is not only 
shifted in frequency but the observed flux density (F Q ) 
is significantly down-boosted compared to its rest-frame 
value (F rf ) with F D = g 3 ~ a x F rf (Lind & Blandford lT3g5)l 
for a power law spectrum of energy index a, which is —1.5 
for the XMM-Newton flare (Porquet et al. l2"UT)3"|l . g 3 ~ a is 
to be averaged over an orbit, indicated by the overline, 
as g varies along the orbit. For a = 0.994, and a disk al- 



most viewed edge-on at an inclination angle of i — 70°, 
g varies between 0.064 and 0.834 at r — r ms and g 3 ~ a = 
0.056, or the radiation received by the observer is reduced 
by a factor of ~18 compared to its rest frame value. For 
larger radii the reduction factor approaches unity, but at 
T = TRmax = 2.5 the radiation is still reduced by a fac- 
tor of 1.7. The reduction becomes much more pronounced 
for disks viewed almost face-on. For i — 15° 0.133 < g < 
0.265 and g 3 ~ a < 0.001. The Lense-Thirring frequency 
and the Keplerian motion at r = r ms have become un- 
observablc like any other frequency close to the marginal 
stable orbit. At r — 2.5, g 3 ~ a < 0.038. Radiation from the 
innermost region is on average dimmed by a factor of some 
hundred, depending on the details of the radial brightness 
distribution. Because of the large reduction factors the 
observation of the inner regions of almost face-on accre- 
tion disks surrounding supermassive extragalactic black 
holes in AGN becomes progressively more difficult with 
increasing a and/or decreasing i. Maybe that this effect 
contributes to some extent to the faintness of the GC black 
hole. 

Although wc have no suggestion which physical pro- 
cess gives rise to the flare we constrain the region from 
where it comes. Since we see the light modulated with fre- 
quencies associated with the marginally stable orbit out to 
the maximum frequency of the radial epicyclic mode and 
in between the bulk of the emission arises from a region of 
1.37 < r < 2.5, i.e. the innermost region of the accretion 
disk. 

4. Conclusion 

We have discovered five distinct periods in the power den- 
sity spectrum of the XMM-Newton X-ray flare of October 
3, 2002. One may wonder about the statistical significance, 
but these five periods, within their measurement uncer- 
tainty, also appear in the power density spectrum of the 
October 26, 2000 X-ray flare observed by Chandra, with 
perhaps one additional period which, however, does not 
show up in the XMM-Newton flare. Further evidence for 
the existence of these periods comes from the June 16, 
2003 infrared flare (Genzel et al. !2003|. The power den- 
sity spectrum shows a clear increase for at least two of 
the periods and possibly for a third period; the remaining 
two periods were not accessible to the infrared observa- 
tions. Each of the periods can be identified with one of 
the four characteristic gravitational modes in accretion 
disks, i.e. Lense-Thirring precession, Kepler motion, ver- 
tical and radial epicyclic oscillation, in such a way that a 
common value for the black hole mass Mbh and the angu- 
lar momentum a emerges, i.e., Mbh = 2.72^°'^ x 10 6 M© 
and a = 0.9939t°;o°74. The available data on M BH derived 
from studies of the orbital motion of the SO-2 star and the 
dark matter concentration in the center of the Milky Way 
are consistent with our result concerning Mbh, which, by 
the way, is independent of the distance to the black hole. 
There are indications of a ^100 s periodicity in the X-ray 
data, both in XMM-Newton and Chandra, but right now, 
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we are reluctant to accept it as firmly established. For 
further progress in the X-ray domain a flare even brighter 
than the XMM-Newton flare is required, and in the in- 
frared band the temporal resolution should be improved 
to ~10s. This is particularly important for simultaneous 
multi-wavelength observations. 
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